The anomalous 14 C-dating f3 decay problem revisited 
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The anomalous inhibition of 14 C-dating f3 decay rate is restudied in terms of shell-model calculations in the jj 
coupling scheme with both realistic and empirical Hamiltonians. It is seen that the accidental cancellation of the 
decay strength is dominated by the mixing effect of two configurations of the final state wave function, \0p~^ 2 ) 

and \0p~} 2 0p~} 2 ) . By decomposing the effective interactions into different tensor components, it is clearly seen 
that the mixing is largely induced by the tensor force. The failure of realistic calculations in reproducing the 
inhibition may be related to its ill description of the monopole component rather than the tensor force. 
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PACS numbers: 21.30.Fe, 21.60.Cs, 27.20.+n 

The anomalously long /3 decay half-life of 14 C has been of 
special theoretical interest since the appearance of the nuclear 
shell model Qly]. The decay involves the J 71 ' = + ground 
state of 14 C and the J n = 1 + ground state of 14 N and satis- 
fies the selection rule for typical allowed Gamow-Teller (GT) 
transitions. However, the extracted transition amplitude from 
experimental half-life is thousands of times smaller than that 
of allowed transitions. The inhibition should be attributed to 
the accidental cancellation of certain components of the in- 
volved state wave functions that contribute to the transition. It 
was recognized that the tensor part of the nuclear force play 
an essential role in inducing the cancellation lUlH]. 

In the original paper of Jancovici and Talmi's yfl, an unrea- 
sonably large tensor force was introduced to induce the can- 
cellation. Later studies show that this problem can be rectified 
by redefining the radial dependence of the tensor component 
(for reviews, see Ref. [2]). One may expect that the shape and 
strength of the tensor force were confined in realistic nucleon- 
nucleon (AW) potentials which are determined by fitting AW 
scattering observables. But the studies of Zamick and collab- 
orators J51-131 showed that the cancellation cannot be repro- 
duced by calculations with realistic Hamiltonians [8] derived 
from microscopic WW potentials like the Hamada- Johnston 
potential and the Bonn potential. This failure was also seen in 
recent calculations of Refs. N9l4llll with modern one-boson- 
exchange and chiral potentials. The problem may indicate 
that the tensor component of the in-medium WW interaction 
is much weaker than that of the bare potential UBl- 

Recently Holt et al. claimed that the problem of realis- 
tic calculations in reproducing the long j3 decay half-life of 
the 14 C can be solved by taking the Brown-Rho scaling in- 
medium modification |9] or three-nucleon corrections II 1 CXI 
of the interaction into account. However, these calculations 
changed dramatically the bulk properties of the effective in- 
teraction and did not shed any light on the role played by the 
tensor force. But perhaps even more serious is that the cal- 
culations fail completely in reproducing the cancellation by 
employing exactly the same approach but with a different WW 
interaction II ill . 

The purpose of this work is to present a general calculation 
on the 14 C j3 decay property with available well-established 



empirical shell-model interactions and compare it with that of 
realistic interactions. By decomposing them into central, spin- 
orbit and tensor parts, we present a quantitative study on the 
role played by different components of effective interactions 
in the cancellation of the 14 C (3 transition amplitude. These 
may help in understanding the problem of realistic interac- 
tions in reproducing the cancellation and in illuminating the 
property the tensor force in the effective interaction and its 
possible in-medium modification effect. It may be also of in- 
terest to mention that the role played by the tensor force in the 
evolution of shell structure has also been extensively probed 
recently (see, e.g., Ref. I121P . 

The wave functions of the 14 C ( J T = + , T = 1) and 14 N 
(J^ = 1 + , T = 0) ground states, denoted by \ipi) and \ipf), 
respectively, can be well described as two holes occupying the 
Opin and 0b 3 / 2 single-particle orbits assuming le O as the in- 
ert core In the jj coupling scheme, the wave functions 
can be written as |5J 



K l°Pl/ 2 2> 



|V/> = a|0p v 2 2 ) + 6|0p3/ 2 0pr/2> + c|0p3/ 2 2 ) , (1) 

where the k and r\ and a, b and c denote the corresponding 
wave function amplitudes. The GT transition matrix element 
is determined by 



M(GT) = %||<7T|hk> 
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+ 2b) + rj(V2b - VZc) 
(2) 

In evaluating above wave function amplitudes and the tran- 
sition amplitude M(GT), one may start from empirical as 
well as realistic shell-model interactions. It has been estab- 
lished that realistic interactions derived from bare WW poten- 
tials are in general close to empirical ones which are obtained 
by fitting experimental data flL3l - [loll . This can also be seen 
from Table U where the diagonal matrix elements of some 
well-defined empirical p-shell interactions Il7l4l9ll and those 
of the realistic interaction of Ref. lUlll are listed for compari- 
son. In fact, the realistic interaction has often been deemed as 
a good starti ng p oint in the construction of empirical Hamil- 
tonians lfl5l Il6ll . It has been shown that, by correcting a 
few monopole (diagonal) interaction terms only, a good de- 
scription on the properties of nuclei over a wide region can 
be obtained by calculations with realistic interactions (see, 
e.g., Ref. lfl4ll for a general description on the property of the 



TABLE I. Comparisons of diagonal matrix elements (ij \ V \ij) JT be- 
tween empirical and realistic interactions. Only the related terms are 
listed for simplicity. The in-medium modified interaction is taken 
from Ref. 1 1 ill . It was calculated with the chiral AW potential assum- 
ing p = po where po is the saturation density. 
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°P3/2 




OPa/2 


OP3/2OP1/2 


OPl/2 


CKQ7] 


-3.19 


-0.26 


-3.58 


-6.22 


-4.15 


HWM Q8J 


-3.68 


-0.15 


-2.62 


-6.55 


-3.95 


WBT Q9J 


-3.85 


-1.22 


-4.16 


-6.86 


-3.45 


WBPQ9J 


-3.91 


-1.15 


-3.86 


-6.94 


-3.45 


Chiral |_1JJ 


-3.28 


-0.61 


-1.19 


-5.20 


-1.67 


Medium Chiral [IT] 


-0.12 


1.91 


1.29 


-3.03 


0.72 



monopole interaction). The origin of this monopole problem 
in realistic NN interactions is not clear yet. It may be related 
to the in-medium or three-nucleon modification of the NN 
interaction which are not taken into account in usual realistic 
shell-model calculations. 

One may expect the in-medium modification was approx- 
imated by simple approaches like the Brown-Rho scaling. 
However, it is known that the Brown-Rho in-medium mod- 
ification may lead to serious changes of the property of the 
effective interaction. The drastic effects of the Brown-Rho 
mass scaling of different baryons on sd-shell effective inter- 
actions were firstly discussed in Ref. ll20ll . It is supported by 
our systematic calculations in different model spaces with var- 
ious AW potentials employing the G-matrix and other renor- 
malization techniques ll2ll 12211 . Calculations of Refs. 191- HTll 
in the p-shell share the same problem, as seen in Table U 
For examples, the interaction terms (0p^ 2 \V\0py 2 ) J=0 ' T=1 

and (0p 2 1/2 \V\0pl /2 ) J=1 > T = change from -0.61 MeV and 
-1.67 MeV to 1.91 MeV and 0.72 MeV, respectively, by ap- 
plying the in-medium modification on the chiral NN poten- 
tial 111 ill . This kind of strange interaction strengths is not sup- 
ported by any existing interactions constructed for this region 
nor experimental data. In fact, it can be easily recognized 
that these matrix elements should be attractive by looking at 
the binding energies and one-neutron separation energies of C 
and N isotopes 12311 . 

The empirical interactions of Refs. lfl7l - [l9ll are constructed 
in the particle-particle channel by assuming 4 He as the in- 
ert core and the single-particle energies as free parameters. 
But in all cases it was assumed that the interactions are are 
same in the hole-hole channel. The energy splitting between 
OPi/2 an d orbits in 15 C and 15 N is calculated to be 

e = e(0p-} 2 ) - e(0p-l) = 6.3 MeV d, 7.3 MeV Hi 
and 6.5 MeV d. The interaction of Ref. Ill is evaluated 
directly in the hole-hole channel with the energy splitting be- 
ing taken as experimental datum, i.e, e = 6.3 MeV. 

The ground state of 14 C is dominated by the configura- 
tion of \0p^ 2 2 ) due to the large spin-orbit splitting between 

orbits 0py 2 and 0p^ 2 . This is supported by calculations 
with both empirical and realistic interactions, as seen from Ta- 



TABLE II. Comparisons of wave functions calculated with empir- 
ical and realistic interactions. All calculations are done with the 
code l24ll except those of lancovici and Talmi's and of the chiral 
potential which are taken from Ref. |3] and Ref. 1 11], respectively. 
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CK|T7J 


0.38 


0.92 


-0.027 


-0.31 


0.95 


HWM (T8J 


0.36 


0.93 


-0.063 


-0.27 


0.96 


WBT [19J 


0.31 


0.95 


0.033 


-0.43 


0.90 


WBP [19] 


0.30 


0.95 


0.014 


-0.41 


0.91 


JT[3J 


0.09 


0.99 


0.20 


-0.41 


0.89 


Zamick [5] 


0.22 


0.98 


0.014 


-0.40 


0.92 


VF[25J 


0.25 


0.97 


0.12 


-0.36 


0.97 


Chiral |JJJ 


0.40 


0.92 


0.14 


-0.68 


0.72 


Medium Chiral [JJJ 


0.26 


0.97 


0.11 


-0.59 


0.80 



ble[TTJ As a result, the coefficient k of Eq. (1) is significantly 
larger than r] (k and r] have the same sign and k = 1 in the 
single-particle limit). The mixing of the two corresponding 
configurations is induced by the non-diagonal matrix element 
(0py 2 \V\0p^ 2 ) J=o ' T=1 for which realistic and empirical in- 
teractions give a similar strength. 

Similarly, one may safely expect that \0p~^ 2 ) should be the 
dominated configuration in the ground state wave function of 
14 N since the other two configurations lie at much higher en- 
ergies. This expectation is supported by all calculations listed 
in Table HD 

Since the amplitude k is much larger than 77, the suppres- 
sion of the Gamow-Teller transition strength [Eq. (2)] should 
be largely due to the cancellation between a and 26 (a and b 
have different signs). That is, the term IQP3/2QP1/2) should 
most likely be the second largest component in the ground 
state wave function of 14 N. Most of our calculations with 
different interactions predict that the absolute value of the 
amplitude of \0p^ 2 ) is small. In this case we should have 

a ~ —(2k + V2i])b in reproducing the cancellation. As 
seen from Table UT1 the problem of the realistic calculations of 
Ref. Ill 111 is that the predicted amplitude a (b) is significantly 
smaller (larger) than expected. 

The ratio between the amplitudes a and b are 
sensitive to the strengths of diagonal interac- 
tion matrix elements (0pf, 2 \V\0pf, 2 ) J=1 > T =° and 

(Op 3 /20pi/2|V r |Op 3 /20p 1 / 2 ) ,/= ^ T=0 and the non-diagonal 
matrix element (0p 3 / 2 0p 1 / 2 |l / |0p^ 2 ) J=1 ' T=0 . If we 
neglect the contribution from the configuration |0p^ 2 ) 
and restrict the calculation to dimension two, the mixing 
between above two components would be solely dom- 
inated by the term (e + (0p 3 / 2 0j?i/ 2 |y |0p 3 / 2 0pi/ 2 ) - 
(0pl /2 \V\0p 2 1/2 ))/(0p 3/2 Pl/2 \V\0p 2 1/2 ). Bearing this in 
mind, one may inspect the in-medium modification calcu- 
lations of Refs. |9 h1 111 . It is interesting to see that the few 
involved non-diagonal matrix elements are not much affected 
by the Brown-Rho or three-nucleon modifications (This may 
be consistent with the fact that the non-diagonal parts of 



TABLE III. The central, spin-orbit (SO) and tensor components of 
the matrix elements {ij\V\kl} of empirical interactions I ljjt figll . 
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1111 


-1.66 


0.58 


0.92 


-1.33 


054 


-0.42 


3333 


-4.43 


0.29 


0.46 


-3.96 


0.33 


-0.21 


1133 


-3.93 


-0.41 


-0.65 


-3.72 


-0.42 


0.29 
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= 1+,T = 






1111 


-4.27 


0.41 


-0.098 


-4.49 


1.11 


-0.075 


1113 


1.08 


-0.21 


0.83 


0.43 


-0.0019 


1.38 


1313 


-5.85 


0.10 


-0.80 


-5.85 


0.35 


-1.36 


1133 


2.72 


-0.13 


-0.46 


1.55 


-0.046 


-0.82 


1333 


3.72 


0.065 


-0.016 


2.41 


0.045 


-0.012 


3333 


-2.96 


0.041 


0.30 


-4.11 


-0.58 


0.53 



modern empirical and realistic interactions are similar to 
each other). For their calculations with the Bonn-B potential, 
which reproduced the cancellation under certain conditions, 
the relative strength of above two diagonal matrix elements 
are rectified at the cost of changing dramatically their absolute 
values 03 El. 

Still it may be of interest to figure out the role played by the 
tensor force in effective interactions in the jj coupling scheme 
(The impossibility of inducing cancellation with effective in- 
teractions without tensor force was first shown by Inglis Oil)- 
Although they are mixed in usual construction of empirical 
interactions, it is possible to separate the central, spin-orbit 
(vector) and tensor force components of the effective inter- 
action through the spin-tensor decomposition procedure l26ll . 
As examples in Table [Til] is listed the different components 
of the effective interactions of Refs. 11181 [19I1 . The decompo- 
sitions of various Cohen-Kurath interactions can be found in 
Ref. 12711 and will not be presented here for simplicity. Similar 
decompositions are also done for realistic interactions derived 
from the state-of-the-art Bonn potential [28]. But as pointed 
out in Ref. l26ll . there is no trivial relation between the spin- 
tensor decomposition of the effective interaction and different 
components of the underlying N N potential. Especially, the 
tensor force of the NN potential may contribute significantly 
to the overall effective interaction when renormalization and 
core-polarization effects are taken into account |8|]. 

In agreement with the earlier calculation of Ref. [27], it is 
seen that the central force components of realistic and em- 
pirical interactions are similar to each other. Noticeable dif- 
ferences are only seen in certain matrix elements of the non- 
central spin-orbit and tensor forces. 

In TablellVlis listed the wave functions calculated with the 
tensor force component removed from the effective interac- 
tions. In this case it is seen that the ground state wave func- 
tion of 14 N is overwhelmingly dominated by the configura- 
tion of \Qp~J 2 ). This is because the non-diagonal matrix el- 
ement (0p3/20pi/2|^|0p^ 2 ) J=1,T=o , which is crucial in in- 
ducing the configuration mixing, is dominated by the contri- 



TABLE IV. Wave functions of C and N calculated with the cen- 
tral force and central and spin-orbit force components of effective 
interactions. 



Interaction 


V 
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Central force only 






CK [IT] 


0.33 


0.94 


-0.15 


-0.085 


0.99 


HWM Q8J 


0.29 


0.96 


-0.14 


-0.086 


0.99 


WBT \19} 


0.30 


0.95 


-0.11 


-0.033 


0.99 



Central plus spin-orbit 



CKQ7J 


0.34 


0.94 


-0.15 


-0.085 


0.99 


HWM Q8J 


0.32 


0.95 


-0.15 


-0.053 


0.99 


WBT [TJ 


0.33 


0.94 


-0.12 


-0.030 


0.99 



bution from the tensor force. As mentioned before, the non- 
diagonal matrix elements are resistant to the Brown-Rho and 
three-nucleon in-medium modifications. Systematic compar- 
isons between realistic and empirical effective interaction in 
this work as well as other studies 014141611 tend to suggest that 
these non-diagonal matrix elements are well defined by the 
underlying NN potential. The problem of realistic interac- 
tions in reproducing the cancellation should be related to its 
uncertainties in the few diagonal monopole terms. 

Summarizing, systematic calculations with a variety of 
shell-model effective interactions predict consistent results on 
the wave functions of the ground states of 14 C and 14 N, both 
of which are dominated by the configuration of 0p^ 2 . It is 

seen that the accidental cancellation of the 14 C-dating j3 de- 
cay amplitude is largely induced by the mixing effect of two 
configurations of the final state wave function, \0p~^ 2 ) and 

| Op^, g QPi/2 ) • T ne mixm g between these two components 
are sensitive to a few diagonal matrix elements and one non- 
diagonal matrix element which is mainly determined by the 
tensor force. The failure of realistic calculations in reproduc- 
ing the inhibition may be related to its ill description of the 
monopole component rather than the tensor force. A rigorous 
correction of the monopole interaction by more comprehen- 
sive in-medium modification or other microscopic approaches 
would be useful. Work in this direction is underway. 
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APPENDIX 

In many cases the wave functions of 14 C and 14 N were cal- 
culated in the LS coupling scheme yfl. The transformation 
between wave functions in LS and jj coupling schemes is 
known in analytic forms in terms of 6j and 9j symbols. To 
facilitate the comparison between wave functions in different 



coupling schemes available on the market, the explicit expres- 
sions for the transformation are listed below as 




_1_ ( 1 y/2 

V3 \ V2 -1 




(3) 



and 
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7! 



/ -1 -4 V10 

v V20 -V5 -V2 
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